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ABSTRACT 


The wind stress calculated from wind velocities measured during 
the Mixed Layer Experiment (MILE) was used as input to a one-dimensional 
wind-driven current model. These model results are compared with ob- 
served currents from the MILE-1 bouy, showing a qualitative agreement. 

MILE was an examination of the upper ocean carried out near Ocean 
weather Station Papa during a 20-day period, August and September of 
1977, which was characterized by two major wind events. 

The observed currents have been analyzed to obtain information about 
their behaviour that could be used in the tuning of the model. For a 
Simulation of the entire Deriod the results are considered only 
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I. INTRODUCTION 


pee eURPOSE OF THE STUDY 

The objective of this research was to investigate how well an one- 
dimensional model of the wind driven current could explain the observed 
near surface current during the Mixed Layer Experiment (MILE). Measure- 
ments of the near-surface currents were taken during MILE, August- 
September 1977, in the vicinity of Ocean Weather Station (OWS) PAPA. 
During this period, the data show the response of the upper ocean to the 
passage of two wind events. 

This study investigates the hypothesis that surface wind can account 
for a large part of the energy and variability of inertial oscillations 
near the ocean surface and that these oscillations are predominantly 
locally generated, as was demonstrated by Pollard and Millard (1970) at 
another site. 

To accomplish this task, a one-dimensional wind-driven current model 
was Solved numerically using a modified "leapfrog" scheme. The modifi- 
cation involved employed current values, u and _ v, which were vertically 
averaged at each time step to a depth equal to that of the observed mixed 
layer, simulating strong boundary layer mixing. The mixing depth used 
was a time-average of the mixed layer depth obtained from buoy temperature 
records and verified with Plessey CTD data taken from the R/V 
OCEANOGRAPHER. 

The data used to check the response of the model were those obtained 


by nine near-surface current meters, at depths from 5 to 32m, of the 





MILE-1 Mooring. The wind stress was computed from wind speed and di- 
rection as recorded by the R/V OCEANOGRAPHER during the experiment. 

It was hoped that the non-wind-driven current would be sufficiently 
small and invariant so that the hypothesis could be tested. 

In the next section, the current meter observations are analyzed by 
using a Fast Fourier Transform. The power svectra of currents from three 
Significant depths illustrate how the currents are dependent on the 
surface wind stress. 

In section three the equations on which the model was based are given. 
It is shown how from those equations a numerical solution is obtained to 
solve the initial value problem that we are dealing with and how the 
model results were computed. 

In the last section the results and conclusions are given. 

If there is no other indication, the abscissa in all figures aives 
time in hours from the start of the experiment, 0500 GMT hours Aucust 


18th. There are 460 hours of recorded currents in all. 


8. IMPORTANCE OF THE PROBLEM 

As Csanady [1981] explains: "The problem of circulation is to de- 
scribe and understand the vattern of the residual or longer-term water 
narticle displacement. The distribution of important water vroperties, 
e.g., temperature, salinity (...}), depends critically on the pattern of 
circulation." The capability of modeling wind-driven currents could 
help to explain one component of that pattern of circulation. 

Modeling wind-driven current is needed in order to improve upper ocean 
thermal structure models, which have great importance in ASW. Fisheries, 
search and rescue operations, and control of pollution in tne ocean are 


other fields where wind-driven currents must be considered. 
8 





C. LITERATURE REVIEW 

Horizontal inertial currents occur in thin layers with a relatively 
reduced horizontal extent, Webster [1968] concludes that inertial 
oscillations are essentially transient phenomena. For a fixed location 
inertial oscillations at different depths may have diverse origins. 
According to Pollard and Millard [1970], Gonella [1971] and Kundu [1976], 
for the surface layer of the ocean, the wind stress is the principal 
mechanism for their generation. 

Inertial currents rotate clockwise (Northern Hemisphere) with a 
period T = n/ (18 sin 6), where @ is the earth's rotation vector and 
8 is the latitude of the observation point. The transmission of the 
movement from the surface to deeper layers may normally be carried out 
either by turbulence in an unbounded homogeneous medium, or by boundary 
effects ina strongly stratified medium. 

To test the dependence of inertial oscillations as a function of 
wind stress the theory of Ekman [1905] was applied in an impulsive 
system. In this theory, the wind provides a stress on the sea surface, 
and the motion is viscously transmitted to lower layers. 

Figures 1 and 2 show four typical examples of the measured vertical 
structures. In figure | both temperature and velocity profiles are con- 
tinuously stratified: the negative gradients of temperature show that 
there 1S no mixing. However, the profiles in figure 2 show a change in 
temperature of less than 0.2° C in the upper 28m, indicating that the 
water is nearly homogeneous and the current is also nearly uniform. The 
more the ocean is stratified, the greater is the tendency for agreement 


with Ekman's model which was developed for a homogeneous ocean. As 
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Gonella [1971] concludes from his work in the Mediterranean this agree- 
ment comes about "because the assumption of a constant coefficient of 
eddy viscosity 1S more realistic in a continuously stratified ocean than 
in an entirely homogeneous ocean where the turbulence may be too intense". 

The approach to the problem in this study is similar to the "slab" 
model of Pollard and Millard [1970] and Garwood [1976], using a depth 
integral of the momentum equation to solve the problem of turbulent 
momentum transfer within the boundary layer. 

The essential principle is that the integrated inertial force, in- 
cluding Coriolis and local time derivative accelerations, just balances 
the applied surface stress less a small damping term. 

To interpret actual current records, two modifications were intro- 
duced in the present work: (i) the momentum equation was integrated to 
a depth equal to that of the mixed layer in each time step, and (ii) a 
linear damping term was applied. By the former procedure, the modeled 
mixed layer becomes unstratified with f = 2n/T as its only natural fre- 
quency. The latter step has the effect of changing that frequency to a 
Slightly smaller value. This may not be a good property when modeling 
the ocean, where inertio-gravity waves have frequencies between f and 
N, the Brunt-Vaisala frequency. Since N is usually greater than f, 
the observed near-inertial currents frequently have a frequency slightly 
greater than f, [Pollard and Millard, 1970]. 

In the one-dimensional model used here surface winds are expected to 
account for a large part of the variability of the energy of the near- 
inertial oscillations near the surface. It was also expected that, as 


in other models [Pollard and Millard, 1970; Kundu, 1976], the changes in 
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phase and amplitude of the forced oscillation could be simulated with 
some success. 

Another problem to be considered is how the energy at inertial fre- 

quencies is removed from the near-surface layers. Pollard [1970], 
Kroll [1975], and Kase [1979] show that one possibility is that energy 
is dispersed through the seasonal thermocline by transferring into in- 
ternal waves in the near-inertial band. Other removal mechanism could 
be the wind [Pollard and Millard, 1970], and in the next section 
evidence will be presented that shows how the wind by, shifting its 
direction, may remove part of the kinetic energy that had been accumu- 

lated earlier. 

Pollard [1980] shows that the amplitude of inertial oscillations can 
fall to zero within a few tens of kilometers, and that the decay of energy 
is due to a slow downward dispersion associated with horizontal spreading, 
Suplemented by the effect of changing winds. In the present work there 
1S particular interest in the downward dispersion and the effect of 


changing wind direction. 





TT. THE MILE DATA 


A. DESCRIPTION 

From 19 August to 6 September 1977, two surface moorings were de- 
ployed in the vecinity of OWS PAPA. The MILE-1 mooring, which provided 
the data used in this study was located at 49° 37'N and 135° 6'W. 
This mooring consisted of an 8m surface toroid anchored with 0.9 scope 
in 4360m of water. Shackled in the mooring line were 19 VACM's of 
which the upper nine were used in this study: those at mean depths of 
pees), 14, 20, 23, 26, 29 and 32m. 

The reason for choosing only those current meters was that the 
mixed layer was at all times more shallow than about 30m, so the measure- 
ments below 32 meter were not considered to be of much significance in 
testing the hypothesis. 

All recorded data were pre-processed and provided on magnetic tape 
in intervals of 15 minutes for a 19-day period, with the following 
exceptions: 

At the 5m level, there are only 5 days of current record. 

At the 8m levle, there are only 16 days of current record. 

At the 14m level there are only 3 days of current record. 

Temperature records at all levels are complete. 

Wind data were recorded hourly aboard the R/V OCEANOGRAPHER. 
Figures 3 and 4 show those wind values. Two strong wind events can be 
seen one strong storm in 22-24 August (starting time 100 hours), and a 
weaker storm in the 31 August-1 September period (starting time 310 


hours). The six-hourly sea-surface pressure maps reveal that both storms 
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were low pressure disturbances which developed rapidly in the Gulf of 
Alaska and moved toward land in an East-Northeast direction. The center 
of the low on 31 August passed directly over station P, causing calm 


conditions for a few hours. That 1s shown by the double pulse in figure 3. 


eee ANALYSIS 

Figure 5 shows hourly temperatures at the studied levels and illus- 
trates the variability of the temperature structure of the layer above the 
seasonal thermocline. Homogeneous layers are indicated in this figure when 
temperature traces from two or more depths coincide. During stratified 
conditions, the vertical or horizontal movement of the water column past 
the sensors causes apparent temperature variability. Most of this vari- 
ability appears to be due to tidal-period internal waves. Above the 
30 meters level, homogeneous and stratified conditions alternate in re- 
sponse to atmospheric conditions, as can be seen in the power spectra in 
figure 5a. 

With the rapid increase of the wind early on 22 August, at about 95 
hours on the time axis, the temperature at 5m dropped and within four 
hours the upper 23m of the water column became vertically homogeneous 
due to strong mixing. During this strong wind event, the mixed layer 
deepened to nearly 35m. After the storm died, by 24 August stratifi- 
cation was reestablished between 17m and the base of the earlier 30m mixed 
layer. During the second storm the weaker wind event of 27-28 August 
mixed the upper layers to about 22m (figure 6). 

Figures 7, 8 and 9 are the power spectra for the u and v_ com- 


ponents of the observed current at three different levels, 11, 23 and 
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32m, for three different time windows. In figure 7 the time window is 
from 0 to 260 hours, only the effects of the first storm are present. 
In figure 9, the time window is from 200 to 460 hours, and only the 
second storm is present. For the spectrum of figure 8, the time window 
is from 100 to 360 hours, and the time period includes both storms. 

Comparison between different depth levels in figures 7 and 9 indi- 
cates different amounts of kinetic energy for each depth at near-inertial 
frequencies. The first storm, as can be seen in figure 4, was charac- 
terized by a pronounced shift from steady easterly winds to strong 
westerly winds. That strong change after a period of at least 100 hours 
of steady winds could be the means by which part of the energy transfer 
between the ocean and the atmosphere could have eliminated the existing 
inertial current before generation of a new inertial current. 

Since the upper ocean 1s assumed to be one-dimensional there must 
be some unsteady mechanism by which the vertical transfer of momentum 
is modulated in time, causing the observed phase and magnitude differences 
between the 32 meter level current and the surface. 

The response of the current in the upper layers to the second storm 
was quite different from the response to the first storm. In this latter 
case the wind shifted from having a northerly component to a strong 
southerly one and, in a few hours, to a northerly component again 
(figure 4). The upper level currents initially appear to be attenuated by 
the southern component winds but then are reinforced when the wind di- 
rection shifts toward the south. Also these winds have a smaller inten- 
Sity and duration, and as it 1S shown by figure 6 the mixed layer did not 


deepen as much as in the first storm. 
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Because the mixed layer had undergone a shallowing process as dis- 
cussed earlier, the transfer of energy to the 32m level was hindered by 
stratification. 

Figures 10 and 11 are the power spectra for the u and v_ compon- 
ents of the current for the total period of observation. The spectra 
for 32m level exhibit a double-peaked feature at the near-inertial period. 
As Pollard [1970] suggests, at that depth it is unlikely that the ampli- 
tude of inertial oscillations can be generated by the action of a single 
storm. 

The time series of observed currents are shown in figures 12 to 19. 
In those figures, the plotted running mean, an average over one inertial 
period, shows the non-inertial current at each level. It should be 
recognized here, that an undetermined part of this non-inertial current 
may be quasigeostrophic and not directly related to the local wind 


conditions. 
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Figure 16. Observed current: 100 to 360 hrs. (u-component). 
Running mean (*), average over 1 inert. period. 
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Figure 17. 
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Figure 18. Observed current: 200 to 460 hrs. (u-component). 
Running mean (*), average over 1 inert. period. 
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Pee DeoCh lr lON Or te MODEL 


A, THEORETICAL EQUATIONS 
In this study, one-dimensional conservation of momentum in the mixed 
layer 1S assumed. This conservation of momentum and the condition of 


incompressibility are reflected by the Navier-Stokes equations of motion, 


Dv Seal 2+ 
ot 2a xve=- 5 Vipee GS cia 
Su. 
ian © 


| 


SeneinG: U = UU. + U + U' 


where u and u' are the horizontal mean and fluctuating components of 
x-component of non-geostrophic velocity, and Ua is its geostrophic 
component; similarly v, v' and ie are the corresponding y-components. 


With the Boussineq approximation and Reynolds averaging, 
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Assuming 
(a) Reynolds stresses are much greater than viscous stresses. 
(b) Vertical velocity field vanishes, w=0. 


(c) Geostrophic flow is separable from wind driven flow, 


and neglecting the horizontal density gradient 
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Then the equations for wind-driven flow become 





Qo 
ct 
i 
h 
<= 
4 
Q 
N 
ag 
= 
a 
— 
= 


Q?> 


Q 
ct] <! 
\ 
tT) 
cS 
\ 
< 
= 
(ae 
Re) 
LJ 


B. NUMERICAL METHOD OF SOLUTION 


To solve equations [1] and [2] a vertical grid from the surface to 


35m, was used. The values of u(z,t) and v(z,t) were then solved 


numerically as an initial-value problem. 


3/ 





Taking the equations for wind-driven flow and applying an eddy 


viscosity closure, 








war. jp alll 
-u'w' = k 77 
Veen av 
-V'W k a7 
and assuming k constant, gives: 
- ne 
au fv + k o4 bey 
OZ 
= oa 
See fie k SS [4] 
dZ 


The boundary conditions are: 


a. At the surface: 





Tar = x HO). AY 
0 OZ Q 
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where T(t) and ty (t) are wind stress in x and y 
directions respectively at time t and o is the density 
of sea water. 


b. At z= -D< -h, a "Slip" condition is prescribed, 


u'w'(-D) 
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viw'(-D) = 0 
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With the eddy viscosity closure, 


3U _ av z 
ay (ab) 2 ne (Ose 


The numerical model used to simulate the current was based on a 
"leapfrog" scheme, with m indicating space (depth) and n_ time, then 
u(z,t) = u(mAz,nAt) = u(m,n). In this manner the component of equations 


[3] and [4] can be written: 


du _ u(m,nt]l) - u(m,n-1) 


e 2At 


ov — v(m,nt1) - v(m,n-1) 








at 2At 
fave FV (iam) 
fu Ss aie 
3¢ k 
kee 5 (u(m-1,n-1) - 2u(m,n-1) + u(m+1,n-1)) 
92Z (Az) 
Nat k 
k —> = (v(m-1,n-1) - 2v(m,n-1) + v(m+tl1,n-1)) 
2 2 
OZ (Az) 


Solving for the (n+1) values, 
u(m,n+1)=Ru(m-1,n-1)+(1-2R)u(m,n-1)+Ru(m+] ,n-1)+2fAtv(myn) 


v(m,n+1 )=Rv(m-1,n-1)+(1-2R)v(m,n-1)+Rv(m+1 ,n-1)-2fAtu(myn) 


where: R = Melee s) ONE. 


At the surface m=M and application of the boundary conditions 
gives: 
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u(M,nt+1) = 2Ru(M-1,n-1) + (1-2R)u(M,n-1) 


+ 2fAtv(M,n) + 2RAz tx(n)/(ko) 


v(M,n+1) = 2Rv(M-1,n-1) + (1-2R)v(M,n-1) 
mee CU eee AZ Ty tly (KO) 


The numerical model is stable when the value of R is equal to or 
less than 0.5. Assuming an eddy viscosity coefficient k of 23 cm/s 
and an increment of depth of Az = Im, and a time step At = 100 
seconds gives an R equal to 0.46. 

The initial conditions were established as those observed values of 
the current speed for each component, u and v, at time 0500 GMT 
August 19th, at the first recorded peak of the v-component of the current. 

The surface stresses TX and Ty were computed every hour as a 


function of the hourly wind speed data for each component Hh and w., 


y 
using 
Wi 
B iE Z 
ae W (wy + Wy) Che 
We 
a 2 Z 
By = Wy (Wy + Wy) Cy? 
with a value for the drag coefficient C, of 1.3 x re. 


D 
Values for the u and v components were calculated for each 


level, from the surface to 35m at intervals of Im, at time steps of 
100 seconds. 

Every hour a new input of wind stress was applied, setting as 
initial conditions for the next run an integrated value for depth in the 


mixed layer. 
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Mixed layer depth as a function of time was deduced from the Mile-1 
temperature time series and checked with Plessey CTD data taken from the 
R/V OCEANOGRAPHER. 

The model was initially run without any damping term, and the 
results were acceptable until the onset of the second storm. At that 
time a clear shift in phase between the recorded and modeled currents 
was observed. In an attempt to avoid a steadily worsening phase problem 
a linear damping term was introduced. According to Pollard and Millard 
[1970], this decay factor models the dispersion effect. Then the 


equations for wind driven flow become: 


ou = a r 1 a 
aa Py a> AU w') = Cu 
dV _ Ge 2) Tae ey 
5t = -fu - ya (Vv W ) Gy 


A damping coefficient of c = 5.03 x ie sec! was chosen from 


Pollard [1970]. Thsi gives an e-folding decay time of 2.3 days. 
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Figure 20 through 22 show the modeled currents compared with the 
recorded ones at three differents levels. 

As was discussed earlier, the dominant frequencies of the model 
and the ocean are slightly different, see figures 23 and 24. There- 
fore when the wind decreases and the current is almost totally inertial 
under small forcing the modeled and observed oscillations will drift 
out of phase. This 1s most apparent in the period before the second 
storm, when the wind speed dropped to a minimum. 

The model does not have the capability to reproduce the observed 
current driven by the sharp change in directions of tne wind during the 
second storm. The modeled and observed currents become out of phase 
at about time 330 hours. After the wind again becomes more steady in 
direction, the modeled and observed currents again agree in pnase. 

Other discrepancies between observed and modeled current could be 
caused either by observational errors or by failure of the model to 
properly treat the physics of the problem, such as neglecting advection. 

It 1s clear that the dissipation of mean kinetic energy is too much 
in the model using R = 1/2.3 days. Figure 22 shows that at the 32m 
level the modeled current is stronger than the recorded one. 

At the other levels the opposite happens. The mode] neecs to be 
tuned to get a better amplitude comparison in the mixed layer. 

Doing this tuning by only comparison with observed currents is 


difficult because of the probable errors in the measurements. As 
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Figure 20. Comparison between observed and modeled currents. 
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Figure 21. Comparison between observed and modeled currents. 
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Figure 23. Power spectrum of modeled current (u-component). 
(Abscissa: c.p.h.; ordinate: (cm/s)°). 
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Figure 24. Power spectrum of modeled current (v-component). 
(Abscissa: c.p.h.; ordinate: (em/s)°). 
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suggested by Halpern et al 11981], there is spurious rotation of the 
rotor of the current meters produced by motions of the surface- 
following buoy mooring. 

Tuning could be done with a different treatment of the eddy vis- 
cosity as well as alternate parameterization of the mean Kinetic energy. 
The first approach would be possible when turbulent viscosity of the 
water as a function of depth and stratification is better known, per- 
mitting detailed quantitative calculations of the vertical stress 
profile. 

Two other factors that could have a relative important effect on the 
currents and are not included due to the limited scope of this study is 
the process of inertial-gravitational wave propagation and the resultant 


horizontal dispersion of mean kinetic energy. 
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